Self-consistent Coulomb picture of an electron-electron bilayer system 
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In this work we implement the self-consistent Thomas-Fermi approach and a local conductivity 
model to an electron-electron bilayer system. The presence of an incompressible strip, originating 
from screening calculations at the top (or bottom) layer is considered as a source of an external 
potential fluctuation to the bottom (or top) layer. This essentially yields modifications to both 
screening properties and the magneto-transport quantities. The effect of the temperature, inter- 
layer distance and density mismatch on the density and the potential fluctuations are investigated. 
It is observed that the existence of the incompressible strips plays an important role simply due to 
their poor screening properties on both screening and the magneto-resistance(MR) properties. Here 
we also report and interpret the observed MR Hysteresis within our model. 

PACS numbers: 73.20.-r, 73.50.Jt, 71.70.Di 



I. INTRODUCTION 

The successful experimental realization of a two di- 
mensional electron system (2DES) has revealed a novel 
technique to exploit the quantum mechanical properties 
of wide range of mesoscopic systems, including integer 
quantized Hal! 1 (IQH) and drag 2 effect. An interesting 
composite two-dimensional (2D) charge system to study 
screening and magneto-transport is the so called bilayer 
system. The basic idea is to bring two 2DES into a close 
proximity, in parallel to each other, perpendicular to the 
growth direction. For such a system it was predicted, 
that transport in the active layer will drive the passive 
layer out of equilibrium. Even if the barrier separat- 
ing the two layers is high and wide enough to prevent 
tunnelling, the inter-layer interactions can still be suffi- 
ciently strong. This effect is known as the drag effect. 
With the improvement of the experimental techniques 
an additional electron or hole layer was also accessible to 
measure magneto-transport quantities of the bilayer sys- 
tems. Motivated by the drag effect experiment a 3 i 4 i 5 i 6 i 7 i 8 , 
both the electrostatic and transport properties of such bi- 
layer systems were investigated theoretically2ii2*iiii2ii&i4 
within the independent electron picture, however, the 
self-consistent treatment of screening was left unresolved. 
On the other hand the direct measurements of the com- 
pressibility of the 2DESi£ had revealed regimes of neg- 
ative thermodynamic compressibility of the interacting 
2DES and were in qualitative agreement with the existing 
theoretical predictions^. Recently, a magneto-resistance 
hysteresis has been reported for the two-dimensional car- 
rier systemsiii^. For the GaAs hole bilayer system, 
Tutuc et al observed hysteretic longitudinal resistance 
at the magnetic field positions where either the major- 
ity (higher density) or minority (lower density) layer 
is at Landau level filling one. They have argued that 
the hysteresis is due to a layer-charge instability, which 
creates domains with different layer densities, whereas 
for the electron bilayer system Pan et al. concluded 
that, the observed hysteresis is due to a spontaneous 



charge transfer via the ohmic contacts. The vast suc- 
cess of the incompressible strip " edge" picture explain- 
ing the IQHE rely on the self-consistent treatment of the 
electron-electron interaction, i.e. screening 19 20 .. One can 
reproduce experimental results of the high precision2i 
quantized Hall (QH) plateaus and the resistance val- 
ues in between these plateaus within a relatively sim- 
ple Thomas-Fermi approximation — and a local version of 
Ohm's lawSi 2 ^ 2 ^. As a result of self-consistent screening 
calculations the 2DES splits into two " domains" , namely 
the quasi-metallic compressible and quasi-insulating in- 
compressible regions. The electron distribution within 
the Hall bar depends on the "pinning" of the Fermi 
level to highly degenerate Landau levels. As expected, 
if the Fermi level is equal to (within few ksT, ks is 
the Boltzmann constant and T is temperature) a Lan- 
dau level with high density of states (DOS) the elec- 
tron system is known to be compressible (locally), other 
wise incompressible. In a very recent experiment sim- 
ilar magneto-transport hysteresis was observed 2 ^ at an 
electron-electron bilayer system. It was attributed to 
thermodynamical non-equilibrium caused by the incom- 
pressible strips similar to the explanation of the hystere- 
sis also observed at single layer systems 2 ^. The exper- 
imental findings were interpreted as the fingerprints of 
dissipationless eddy currents driven by an induced gradi- 
ent in the electrochemical potential over the incompress- 
ible regions. Huels et al. concluded that, when measuring 
the Hall resistance while sweeping the magnetic field, the 
2DES cannot be considered as being at thermodynamic 
equilibrium at the plateau regimes, where one expects 
incompressible regions. 

In this paper we apply the self-consistent scheme de- 
veloped in the previous works 2 ^ 2 ^* 2 ^ 2 ^ to an electron- 
electron bilayer system. For this model system, we can 
investigate both self-consistent screening and magneto- 
transport properties within the linear response regime 
under QH conditions. We show that the existence of the 
incompressible regions in one of the layers, affects the 
other layer density profile strongly by creating poten- 
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tial fluctuations. We observe that these potential fluc- 
tuations modify the magneto-transport quantities. Here 
we present our model in detail, that explains some of 
the recent experimental findings— , where the magneto- 
transport hysteresis was observed for mismatched den- 
sities, whenever one of the layers is in the QH plateau 
regime. 

The organization of this work is as follows: first the 
bilayer geometry and fixed external potentials acting on 
the electron layers are introduced. Then we discuss the 
electron density and electrostatic potential profiles of 
the system, which are obtained within the self-consistent 
Thomas-Fermi-Poisson approach--. Second, we system- 
atically investigate the influence of the incompressible 
strips depending on temperature, inter-layer distance and 
the density mismatch, then shortly discuss the effect of 
density mismatch on both electrostatic and transport 
properties. Here, we use the scheme presented by Sid- 
diki and Gerhardts (SG)24, to calculate the Hall and the 
longitudinal resistances. We examine the effects of the in- 
compressible strips from the point of external potential 
fluctuations based on the arguments given in the same 
work. There they have concluded that a perturbing ex- 
ternal potential may shift, widen and/or stabilize the QH 
plateaus. 



II. 



THE MODEL 



In a typical electron bilayer sample, a silicon doped 
thick (AlGa) As layer is grown on top of a GaAs substrate. 
This is followed by the bottom GaAs quantum well, 
separated from the top 2DES by an un-doped AlGaAs 
spacer. On top of the upper 2DES again a silicon doped 
(AlGa)As layer is grown, which is capped by a GaAs 
layer. So that two 2DES are placed in a close proxim- 
ity, which are confined by remote donors and the sample 
is capped by top and/or bottom gates that controls the 
electron density of each layer. In related experiments 
both gates are used to tune the electron densities from 
top and bottom layers, whereas the gate potential profile 
depends on the sample geometry and the applied gate 
bias. The electrons are symmetrically (with respect to 
the growth direction) confined by AlGaAs layers each of 
which contains a plain with Si doping (5— doping), and 
are separated by a spacer of thickness (h ~ 10 — 30nm). 
For such a separation thickness the bilayer system is 
then known to be electronically decoupled and in non- 
equilibrium, i.e. can be represented by two different elec- 
trochemical potentials. Electron tunnelling between the 
layers is not possible. 

In Fig.^a schematic drawing of a mesa etched bilayer 
system is shown, which consists of two 2DESs (minus 
signs) , two donor layers (plus signs) and a top gate (gray 
area). We model the bilayer system such that the bot- 
tom 2DES lies on the z = plane with a number den- 
sity nfi(x) in the stripe — d < x < d, where 2d is the 
width of the sample and the top layer is in the plane 



z = h, with the electron density n^(x). Relevant to the 
considered experiments, we assume that the donors (or 
equivalently the background charges) are symmetrically 
separated from the electron layers placed at z = — c and 
z = c + h, for bottom and top layers respectively, having 



a constant surface density n. 



uq. Here we also 



assume translation invariance in the y direction. The 
electron density of the top 2DES is governed by the top 
gate, located atz~c + h + f = z g and the electron 
channels are formed in the interval Ixl < b. One obtains 
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FIG. 1: The cross-section of the bilayer geometry. Sample is 
capped by a gate from top in order to change the electron den- 
sity of the top-layer. Top (bottom) 2DES lies below (above) 
the top (bottom) layer donors at a distance, c. Electron layers 
are separated from each other by a dielectric spacer having 
thickness, h. 

the total electrostatic potential of an electron on the line 
(x, z) due to a line-charge at (xq, zq) from 22 



V(x, z) = Vb g (x, z) + V H (x, z), 



(1) 



and 

2e 2 f Xr 

V H (x,z) = — dx K(x,xo 7 z,z )n ei (x ,z ) 7 (2) 
K Jx, 

where — e is the charge of an electron, R an average back- 
ground dielectric constant, and the kernel K(x, xq, z, z ) 
solves Poisson's equation under the relevant boundary 
conditions given by 



K(x,x ,z,z ) 



COS 2 £,(x + X ) + 7 2 



"2 ln L - ,■„) + -.* 



(3) 



where the z— dependence is given by 7 = sinh(7r|z — 
Zo|/4d), due to a missprint a factor 1/2 was missing in 
our previous work, although the numerics included this 
factor. The confining (background) potential is obtained 
by inserting a constant number density (no) of the back- 
ground charges into Eq. @. It is assumed that the gate 
can be described by an induced charge distribution n g (x) , 
residing on the plane z g . Then, from Eq.© the gate po- 
tential can be written as 



2e 2 f d 

V g {x,z) = —r- I dxoK(x,xo,z,z g )nJxo). 
« J-d 



(4) 
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In order to obtain a flat (gate) potential profile at the 
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FIG. 2: The gate potential profile (solid-line) at z = ft, 
together with the generating positive charge distribution 
(dotted-line) against position. The additional charges resid- 
ing on the gate lead to a stronger confining potential in the top 
layer, resulting in a higher average electron density. The gate 
potential strength is taken to be n®/no = 0.02, where Eq is the 
pinch-off energy. It is defined as Eq = (2-Ke 2 nod/R)(ii 2 /diG) 
normalized with G = 0.915965594, the Catalan's constant— 
for 7 = 0. The donor density is no = 4 • 10 11 cm -2 . For other 
parameters see text. 



bulk we choose the induced charge distribution as 
n g (x) = n° g (l + a(x/d) 2 ), 



(5) 



where n° determines the strength of the gate potential 
whereas a gives the slope of the induced charge distri- 
bution. In Fig. |21 the distribution of the positive charges 
and the resulting potential profile is shown for a = 0.7. 
One can add more electrons to the top layer by setting 
rig/no to a positive value, while keeping the depletion 
length d—b fixed at zero temperature and for vanishing 
magnetic field. With such treatment of the gate, the av- 
erage electron densities are changed while the depletion 
length is kept constant. Similarly the Hartree potential 
is calculated from Eq.@ for the top layer as 



2e 2 



-d 



dx K(x, xq, z, z a = h)n^i(x ) 



and for the bottom layer 

2e 2 f d 
V£(x,z) = — / dx K(x, 
« J-d 

Then Eq.Q can be rewritten as 



x ,z,z = 0)7i2(x ). 



V(x, z) = -V b T g (x, z) - V b B g (x, z) - V g {x, z) 



+V£(x,z) + V*(x,z). 



In the next step the electron densities are calculated 
within the TFA: 



T,B, 



'(x) = f dED(E)f([E + V(x,z T , B ) - f4,B\/k B T), 

(7) 

with D(E) the Landau density of states (DOS), f(E) = 
X/[exp(E/kBT)+l] the Fermi function, fij, B is the chem- 
ical potential (being constant in the equilibrium state) 
and zt.b the position of top and bottom layer, respec- 
tively. We will use the Landau DOS described by 



1 °° 

D(E) = ^Y, 6 ( E - E ^> 



(8) 



n=0 



as default, unless other definitions are given. In our nu- 
merical calculations, we start with zero temperature and 
magnetic field and initially obtain the confining potential 
created by its own donors for each layer. Then, we obtain 



LU 



-0.25 



-0.5 



-0.75 



-1 



\\ 
\\ 


- V bg T (x,h)=V bg B (x,0) 


i\ 


\ \ 

\ A 
\ \ 
\ \ 


- V bg B (x,h)=V bg T (x,0) 


// / 

ft 
// / 

// 


- \ "\ 
\ *\ 

\ Vv 

\ \\ 

\ V\ 


— Oo 


// 
// I 
// / 
// 1 
// I 

V / 









-1500 



-500 



500 



1500 



X fnm) 



FIG. 3: Total background potential created by the top and 
bottom layer donors (thick solid-line). The potential contri- 
butions to the bottom 2DES (z = 0) come from top donor 
layer (V b g(x, 0), thin solid line) and from bottom donor layer 
(Vt, g (x, 0), dashed line). Numerical calculations are performed 
for GaAs heterostructure with the parameters h = 15nm and 
c = 60nm. The numbers of donors are fixed to 4 x 10 11 cm -2 



per layer, for a sample width 2d 
pinch-off energy Eq. 



3/xnm, with an average 



the electron densities and the Hartree potentials of both 
layers just as in the case of a single layer—. Knowing 
the electron distribution of the top (bottom) layer via 
Eqs.lQ) and J2J we calculate the potential acting on the 
bottom (top) layer similar to the single layer case. So 
that the total electrostatic potential of an electron at the 
top layer is 

V T (x) = V ext (x, h) + V§{x, h) + V£(x, h), (9) 

where the total fixed external potential acting on the 
layers is V ext {x,z) = V b t ° t (x,z) + V g (x,z) and we de- 
(6) note the total background potential by V b t ° t (x,z) = 



4 



Vbg( x ' z ) + v bg ( x > z ) ( see FigEJ. Note that the first two 
terms in Eq. ©, are considered to be external and the 
last term is the intra-layer Hartree potential. Similarly 
for the bottom layer the total electrostatic potential is 

V B (x) = V ext {x,0) + V£(x, 0) + Vg(x,0). (10) 

In each iteration step an accurate numerical conver- 
gency is achieved for a single layer, then the Hartree po- 
tential is added to the other layers external potential. 
Intra-layer self-consistencies are obtained by Newton- 
Raphson iteration and inter-layer self-consistency is ob- 
tained by direct iteration. One can then use this so- 
lution as an initial value and obtain the density and 
screened potential for finite magnetic field and tempera- 
ture. We scale energies by the average pinch-off energy 
i?o = {Eq + E§)/2 = E , since we always assume sym- 
metric donor distribution, e.g. £l c — huj c /Eo (and set the 
gate voltage to be zero, thus no extra charges reside on 
the top gate). The lengths are scaled by the screening 
length a (= a* B /2) expressed in terms of effective Bohr 
radius, a* B = Rh 2 /(me 2 ). The electron density and the 
electrostatic potential can be calculated self-consistently 
by the above scheme within the TFA. 



III. SCREENING RESULTS 

In this section we present our results calculated within 
the self-consistent scheme, starting by discussing the zero 
temperature, zero magnetic field limit and then compare 
the obtained density profiles with the finite temperature 
and field profiles. The aim of this investigation is to clar- 
ify the effect of the quantizing perpendicular magnetic 
field, which introduces local charge imbalances due to 
formation of the incompressible strips. We have already 
seen that the electron density distribution is highly sen- 
sitive to the applied external magnetic and electric fields. 
Therefore even very small changes in these external pa- 
rameters affect the density and potential profiles drasti- 
cally. 

It is well known 22 that the formation of the compress- 
ible and incompressible strips results in an inhomoge- 
neous density distribution that deviates from the zero 
field profile. This deviation creates a local charge imbal- 
ance generating a potential fluctuation. We describe this 
fluctuation by, 

AV(x) = [V(x, 0, 0) - V(x, T, B)]/Sh, (11) 

where the self-consistent potentials are calculated at zero 
and finite magnetic field and temperature, respectively 
(see FigEb). Here after Q 2 (= hu c /E = 0.328 x 10~ 2 ) 
represents the dimensionless cyclotron energy at aver- 
age filling factor two, since we always consider situa- 
tions, where average filling factor is around two. We 
start our analysis by discussing the effect of the lo- 
cal charge deviation from its equilibrium distribution at 



T = , B = 0. Figure UJi depicts the electron densi- 
ties calculated within the SCTFPA for vanishing and fi- 
nite magnetic field and temperature, where we used 300 
mesh points to span a single layer. The typical parame- 
ters used in the calculations are: d — 1570nm, c = 60nm, 
h = 15nm, / = lOOnm, hq — 4 x 10 n c:m -2 and for un- 
biased gate = = 3.31 x 10 11 cm -2 which result 
in a Fermi energy (Ep) ~ 14meV. A positive (with re- 
spect to the electrons) potential bias is applied to the 
top gate [for details check Eq.© and related text] so 
that more electrons are populated to the top layer re- 
sulting in a density mismatch. The curves for finite field 
and temperature show a considerable deviation from the 
curves for zero field and temperature in the intervals 
(200nm< \x\ < lOOOnm), where one observes incompress- 
ible strips at the top layer. In the inset we concentrate 
on this interval. In Fig. 0t we compare the T = 0, B = 
density curves to the T ^ 0, B ^ ones. In the inter- 
val — 600nm< x < — 300nm there are less electrons at the 
top layer due to the formation of an incompressible strip. 
This yields a less repulsive inter-layer Coulomb interac- 
tion. So that more electrons are populated locally at the 
bottom layer. Similar arguments hold for the left-hand 
side of the incompressible strip (— 800nm< x < — 600nm) 
but now the potential becomes more repulsive thus more 
electrons are depleted from the bottom layer locally. The 
corresponding potential variations are shown in Fig^D, 
the peak like behavior near the edges results from tem- 
perature difference, whereas the structure observed at 
the top layer is due to the formation of the incompress- 
ible strip in the same layer. Since the bottom layer is 
completely compressible the potential variation does not 
show any non-monotonous feature and the fluctuation is 
perfectly screened. A relatively large variation is seen 
within the depletion regions. The reason for this is that, 
at finite temperature the density profiles leak out at the 
edges. At very strong magnetic fields (is(0) < 2) only the 
lowest Landau level is partially occupied, i.e. high DOS, 
thus the density profiles for vanishing and for very strong 
magnetic field should look very similar. However, for in- 
termediate B there is a difference due to the quantiz- 
ing magnetic field, that creates dipolar (incompressible) 
strips at the top layer, which have an influence on the 
bottom layer via the strong Coulomb interaction. There- 
fore, it is essential to examine the formation of the incom- 
pressible strips at strong magnetic fields. This is done by 
manipulating the sheet electron densities by applying a 
finite gate bias. The gate controls the existence and po- 
sitions of the incompressible strips indirectly, so that one 
can examine the screening effects of the incompressible 
strips considering the electron density distributions. In 
the later sections, potential fluctuations created by these 
local charge imbalances, i.e. the incompressible strips, 
will be connected to the magneto-transport quantities, 
where we explicitly show the impact of the incompress- 
ible strips on the Hall resistances. 
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FIG. 4: [a] Electron densities for a finite magnetic field (such 
that, z/*°* = 1.6) and at default temperature (kT/Eo = 5 x 
1(T 5 or kT/E F ~ 0.01), for top and bottom (thin solid-line) 
layers. Also for vanishing field and temperature (dotted lines), 
[b] The resulting potential variation due to the formation of 
the incompressible strips, where the superscripts refers to top 
(T) and bottom (B). The depletion length is set to be 150nm 
and the density mismatch is governed by applying a finite gat 
potential Vq/E^ = 0.02, resulting in ftj = 3.42 x 10 n cm" 
and = 3.18 x 10 11 cm~ 2 . The inset shows the region whei 
an incompressible strip exists at the top layer. 



A. Intra-layer distance, temperature and density 
mismatch 

It is well known that the mutual Coulomb interat 
tion is a strong long-range interaction. Thus, a chang 
in the charge distribution, compared to the equilibriur 
distribution (at T = 0,B = 0), produces a considei 
able effect on the observable quantities even at large dis 
tances. In the previous section it is shown that such 
local charge imbalance, connected with a potential flue 
tuation, is created due to formation of the incompresf 
ible strips. Here we investigate the density distribution 
of the layers by applying a positive gate voltage, henc 
populating the top 2DES and vary the intra-layer dis 
tance and temperature. The average total filling factoi 
^av = i vT + is kept constant and the evolution 

of the incompressible strip at the top layer, and its effect 
to the bottom layer, is examined. In Fig. [S] we show the 
local filling factors of both layers versus position, where 
the top layer (upper set of curves) has an incompressible 
strip in the interval — 850nm< x < — 600nm for different 
inter-layer distances. The influence of the incompressible 
strip in the top layer on the electron distribution of the 
bottom layer disappears rapidly although the intra-layer 
distance is changed rather smoothly. This is due to expo- 
nential decay of the amplitude of the Coulomb potential 
in the z— direction, i.e. V(q,z) ~ V^exp— q\z\. The 
effective confining potential (ECP) experienced by each 
layer, of course, depends strongly on the inter-layer dis- 
tance, hence for each h value the number of electrons at 
each layer changes. This is depicted at the inset of Fig0 
there we concentrate to the center filling factor center. 
We observe that by increasing h(= 30,60 and 120nm) 
the number of electrons at the center decreases (at fixed 
magnetic field) and the ECP becomes less confining thus 
a flatter density distribution is observed. Interestingly, 
for the largest separation (h = 200nm) the bottom layer 



becomes widely incompressible at the bulk, hence semi- 
transparent, and the electrons of the top layer start to 
see the donors of the bottom layer. This specific config- 
uration leads accumulation of the top layer electrons to 
the bulk and is a clear indication of non-linear screening, 
due to incompressible strips, in such a bilayer system. 
This fact emphasizes that even small a change at h has a 
strong influence on the electron density profile. The ex- 
ponential decay of the Coulomb interaction in the growth 
direction essentially determines whether the two 2DESs 
are strongly coupled or not, depending on the formation 
of the incompressible strips, i.e. screening. In the rele- 
vant experiments the distance between the two layers is 
fixed during the growth process. Hence we proceed our 
investigation by fixing h = 15nm and vary the electron 
temoerature. For such a seDaration the svstem is known 




X(nm) 

FIG. 5: Local filling factors versus position at a finite mag- 
netic field and for h = 30, 60, 120 and 200 nm, at default 
temperature. The top layer (upper set) exhibits an incom- 
pressible strip and whereas bottom layer (lower set) is com- 
pressible all over the electron channel, except h = 200nm. 
The inset depicts the center filling factor of the top layer for 
four h values. 

In Fig. [S] the temperature dependence of the incom- 
pressible strip residing in the top layer and the com- 
pressible region at the bottom layer is shown. At sin- 
gle layer geometries it is well knownitLSi that the quan- 
tizing effect of the magnetic field becomes inefficient, if 
the thermal energy of the system becomes larger than 
few percents of the cyclotron energy. Here we also ob- 
serve a similar behavior for the bilayer system where 
the incompressible strip at the top layer and, moreover, 
the local charge imbalance seen at the bottom layer dis- 
appears while increasing the dimensionless temperature 
t{= kT/Ep) from 0.01 to 0.08. This fact shows that the 
local inhomogeneity at the bottom layer density distri- 
bution is only due to the dipolar strip at the top layer 
and is very sensitive to the temperature and exists for 
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t < 0.04. 
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FIG. 6: Density profiles across the incompressible strip at 
the top layer as a function of temperature. For a fixed density 



mismatch, 



0.93, at v, 



tot 
av 



1.6 with h — 15nm. 



Relevant to the experiments considered, density mis- 
match is another tunable parameter, like the tempera- 
ture. For this purpose now we fix the temperature and 
intra-layer distance and vary the potential of the top gate 
essentially by changing the number of induced positive 
charges. Figure presents position dependent filling fac- 
tors for such a positively charged gate, simulating differ- 
ent density mismatches by applying a gate bias voltage 
Vq. For a slight mismatch one does not observe a promi- 
nent change, except more electrons are accumulated at 
the center of the top layer. Due to the stronger Coulomb 
repulsion the bulk of the bottom layer is more depopu- 
lated and shows a flatter profile ( see FigEJi). In FigQj, 
this feature is more pronounced at a higher gate bias and, 
in turn, the bottom layer is forced to be incompressible 
at the bulk, leading the effective external potential to 
be more confining for the top layer. The outer edge re- 
gions of incompressible strips residing at the top layer 
suppress the electrons beneath, at the bottom layer. In- 
creasing the density mismatch in favor of the top layer, in 
Fig|7Ji and d, it is observed that more electrons start to 
accumulate at the bulk of the top layer and small density 
fluctuations can be seen at the bottom layer, however, we 
do not see any prominent change at the density profiles 
since the bottom layer can screen perfectly. We close 
our short discussion by noting that, even a small amount 
of density mismatch can lead to a drastic change in the 
density profiles. Moreover this change is enhanced by 
the existence of the incompressible strip, which results 
from the non linear screening in the presence of exter- 
nal magnetic field. This high sensitivity to the external 
electric field (here the gate) is clearly seen, if the poten- 
tial profiles are considered. We shall emphasize, that a 
wide incompressible strip formed at the bulk of one of the 
layers does not necessarily imply that this layer becomes 
completely transparent to its donors, since the electrons 



residing in the incompressible strip create an electric field 
which can still partially cancel the electric field generated 
by its donors. 
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FIG. 7: The top layer is populated by applying a finite gate 
potential Vo/Eq [a] 0.01, [b] 0.02, [c] 0.03, [d] 0.05 correspond- 
ing to density mismatches of nfi/n^i = 0.96,0.93,0.89,0.84, 
respectively. Given at default temperature, with the depletion 
length, b/d = 0.9 and v l °^ = 1.67. 



B. Potential fluctuations 

In this section we examine the effects of the incom- 
pressible strips in one layer on the external potential 
profile of the other layer by comparing what we call "in- 
teracting" and " non- interacting" systems. In both cases 
we start with the self-consistent calculation of the den- 
sity and potential profiles at T = 0, B = 0. Then we 
focus on one layer, bottom or top, which we call the "ac- 
tive" layer while the other layer is called the "passive" 
one. Now we keep the T = 0, B = density and Hartree 
potential profile of the passive layer fixed and calculate 
in the corresponding, B— independent external poten- 
tial of the active layer its density and potential profile 
self-consistently at finite T and B. This yields the to- 
tal potential Vn(x 7 z) of the non-interacting system. It 
takes into account the B— dependent intra-layer screen- 
ing properties in the active layer, but not the B— depen- 
dent changes of the intra-layer screening in the passive 
layer and of the inter-layer screening. Finally we drop the 
restriction on the passive layer and calculate density and 
potential profiles for both layers at finite T and B fully 
self-consistently. This yields the total potential Vi(x,z) 
of the interacting system. 

The potential variation (at finite temperature and mag- 
netic field, indexed by F as a subscript) 



A F V(x) = [Vr(x, z) - V N {x, z)]/fi 2l 



(12) 
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taken at the z— value of the active layer, describes the 
B— dependent change of the interaction of the active 
with the passive layer. It is suitable for studying the ef- 
fect of the incompressible strips in the passive layer on 
the effective potential in the active layer. Fig. [S] shows 
potential variation and filling factors across the sample 
at four different magnetic field values where the super- 
script T(B) indicates that the top (bottom) layer is the 
active layer. In Fig.|5k the magnetic field strength is cho- 
sen such that both layers are compressible, i.e. the cen- 
ter filling factors of both layers are slightly below two. 
The potential variation shows a characteristic behavior 
as the layers are both compressible all over the sample, 
the screening is nearly perfect and we do not observe 
any significant potential fluctuation. If one decreases the 
magnetic field strength and obtains a wide incompressible 
strip in the bulk of the top layer (Fig. |SJd) , a large poten- 
tial variation is observed at this layer. This is nothing but 
the charge quadrupole at the center, treated by Ref.[3(j 
for a single layer geometry. The variation drastically in- 
creases and becomes almost ~ %30f^2- Meanwhile, the 
variation of the bottom layer does not show any signif- 
icant change. The explanation of these observations is 
twofold; first as an incompressible strip is formed at the 
top layer, the finite field density profile strongly devi- 
ates from the zero field profile, which essentially creates 
a huge local (within the incompressible strip) external po- 
tential fluctuation to the bottom layer. Second, the bot- 
tom layer is completely compressible, so that this large 
potential fluctuation can be screened by redistribution of 
the bottom layer electrons, resulting in a density devia- 
tion from the zero B profile. Note that as a consequence 
of the self-consistency the deviation is spread all over the 
bottom layer, which also generates an external potential 
fluctuation to the top layer. Here, as we examine how the 
top layer responds to this external potential fluctuation, 
we remind the reader that there are two different regions 
with different screening properties: (i) The compressible 
regions (0.5/im< \x\ < 1.5/xm interval in FiglSJ}) show 
similar features to the bottom layer and the external po- 
tential fluctuation is well screened, (ii) The incompress- 
ible region (0/mi< \x\ < 0.5/zm interval in FiglSJ)) at the 
bulk can poorly screen the fluctuation and one observes 
that ApV(x) can be large as half the Fermi energy. In 
Fig. [HJ;, the magnetic field strength is slightly decreased 
in order to obtain a wide incompressible strip at the bot- 
tom layer, meanwhile the wide incompressible strip of 
the top layer splits into two ribbons which are shifted to- 
wards the edges. For this B value, one observes both the 
quadrupole (for the bottom layer) and the dipole (at the 
top layer) moments. Figure [SJi shows the case, where 
two well developed incompressible strips are present at 
both of the layers and the potential fluctuations, created 
by the charge dipoles, are confined to these regions. Dif- 
ferent than the previous case we only observe one side of 
the dipole moment, since the other half is screened by the 
other layer. We should emphasize that the variation am- 
plitude depends on the width of the incompressible strip 




X (urn) X (urn) 



FIG. 8: Position dependent filling factors and magnetic field 
induced potential variations, for P*? [a] 1.44 [b] 1.59 [c] 1.69 
[d] 1.84. For clarity all the potential variations are shifted by 
an amount of 1.5. The gate potential is set to Vo/Eq = 0.02, 
for fixed depletion and at default temperature. 

which generates the fluctuation, since the charge imbal- 
ance becomes larger if the incompressible strip is wide. 
The finding is simply that the potential fluctuations in 
a layer exist only within incompressible regions and are 
screened in compressible ones. 

C. Modulated donor distribution 

It is rather a common technique, to add an external 
(harmonic) modulation potential to the confinement po- 
tential in order to study the peculiar low-temperature 
screening properties of a 2DES, which has been done for 
a single layer geometry22i2ii22. Although a double layer 
geometry is a very promising system, even a qualitative 
investigation of such a modulation is still missing in the 
literature. In this section we provide a simple model to 
discuss the effect of this modulation on screening and 
compare our results, qualitatively, with the single layer 
ones. 

A harmonic density modulation is added to the back- 
ground charge distribution in order to obtain a potential 
modulation. Here we should also note that such a modu- 
lation will be used to simulate the long-range fluctuations 
of the confining potential in the following sections^. It 
is assumed that the spatial distribution of the donors is 
given by 

n m {x) — no(l — i5cos(A7r(a;/2d))) (13) 

where A is an odd integer preserving the boundary condi- 
tions and 5 describes the strength of the modulation. In 
the following discussion we will fix the density mismatch, 
nfl/n^i = 0.77, modulation period, A = 11, temperature, 
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t = 0.01, and the average filling factor, = 1.34, to cal- 
culate the electrostatic quantities. We see in Fig|5]that 
the top layer is in phase with the modulation, whereas 
the bottom layers phase is shifted by an amount of 7r/2, 
where the modulation strength is chosen to be one per- 
cent. For the given parameters both layers are compress- 
ible and screening is still linear (see FigllOlregion (i)), i.e. 
the electrons can redistribute according to the applied 
external potential. A similar case for single layer geome- 
tries has been extensively studied in Ref.:.22l33| and the 
linear dielectric function is given by e(q) = 1 + l/a q 
(~ 19 for our parameters). This linear screening approx- 
imation breaks down when the amplitude of the screened 
potential becomes equal to the Fermi energy of the un- 
modulated system^. From this we would expect that 
the breakdown amplitudes of the two layers should be 
directly proportional to the density mismatch. In order 
to test this and examine the non-linear screening regime 
we increase the strength of the modulation monotonously 
and look at the variation of the screened potential defined 
by, 



var\V 



T.Bi 



V T - B {x = 0) 



V T - B (x 



0.3/xm)]/O 2 . 

(14) 

Figure 110b presents this variation, for top (solid line) 
and bottom (dashed line) layers. In the regime denoted 
by (i) both layers are compressible and the density pro- 
file can be characterized similar to Fig|5] With increas- 
ing the modulation amplitude first the top layer enters 
to the non-linear screening regime, since incompressible 
strips are formed (see Fig llOb h This is pronounced as a 
jump in the variation an saturates at 8 ~ 0.05 (regime 
(ii)), meanwhile the bottom layer is still compressible, 
(i.e in the linear regime) and compensates (screens) the 
potential fluctuation generated [cf. discussion related to 
Fig|Ht>] • From the saturation point we can easily predict 
the similar point of the bottom layer to be S ~ 0.07. Ac- 
tually we see that regime (iii) of FiglTUk starts at the ex- 
pected value, where incompressible strips are present at 
both layers showing a characteristic distribution similar 
to Fig|SJ;. Here screening becomes linear again, but now 
with e(q) = 1 + (Dt / Dq) / (aoq) , where Dt is the thermo- 
dynamic DOS in a Landau Level, which can be estimated 
by D t /Dq ~ huj c /(4k B T). For 5 > 0.15 the bottom layer 
is split to five narrower channels (see FigllOtO and the 
linear screening is broken around S ~ 0.16. From the 
analytical expressions given in Ref . |32| one can estimate 
the breakdown amplitude as 



2E§ 



Jbd 



E 



0.16, 



(15) 



which indeed is in very good agreement with the obtained 
numerical result. The small difference is due to the higher 
temperature and the decay of mutual Coulomb interac- 
tion in growth direction, likewise the difference at the 
slope of the plateaus. We also observe that the plateaus 
occur at the integer multiples of the individual filling fac- 
tors, similar to the single layer system. Since we start 



with a situation where the average filling factors of the 
both layers are below two we see only one plateau, as 
discussed for the single layer geometries the number of 
the plateaus depend on the average filling factor without 
modulation. Finally we would like to note that the linear 
screening will break down for the top layer for S ~ 0.2 
which can seen easily from Ea ll5l 
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FIG. 9: Spatial distribution of the donors (solid line), to- 
gether with the local filling factors of top (dashes) and bottom 
(dashed-dots) layers at default temperature. 
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FIG. 10: [a] The variation of the screened potential versus 
the modulation amplitude, together with the typical electron 
densities corresponding to different screening regimes [b] ii, 
5 = 0.026, [c] iii, 6 = 0.073 and [d] iv, 6 = 0.157. For top 
(solid lines) and bottom (dashed lines) layers. 

The potential fluctuations and the donor modulation 
discussed above plays an important role when one con- 
siders a fixed external current flowing through both of 
the layers. These fluctuations are generated by local 
charge imbalances, with respect to zero field density dis- 
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FIG. 11: Comparison of the interacting (solid-lines) and non- 
interacting (dashed-lines) Hall resistance curves for fixed elec- 
tron densities at default temperature. The Hall resistance 
curve is shifted for the top layer as the potential fluctuation 
created by the bottom layers incompressible strip. 

tributions, and can be observed in the interval where the 
other layer has an incompressible strip as they are poorly 
screened. 



IV. THE HALL RESISTANCE CURVES 

In this section we investigate the effects of potential 
fluctuations on the Hall resistances using the findings of 
within the linear response regime. In their work 
it was concluded that the (long-range) potential fluctu- 
ations can widen, stabilize and shift the quantized Hall 
(QH) plateaus as they affect the position and the exis- 
tence of the incompressible strips^. In addition to the 
long-range fluctuations, in this section we also include the 
potential fluctuations generated by the local charge im- 
balances. Here we use the general expressions derived in 
Ref.[24|, for a given electron density and fixed current, to 
calculate the Hall and longitudinal resistances, of the bi- 
layer system. We also note that the averaging procedure 
of the conductivities is carried over a length scale (A a „), 
which is comparable with the Fermi wavelength and, in 
particular is set to be 40nm in our calculations. We ap- 
ply these results to two cases considered in Sec. IIII Bl 
namely the interacting and non-interacting systems. We 
remind that, in the case of non-interacting system, the 
active layer does not have information about the density 
inhomogeneities caused by the incompressible strips of 
the passive layer. Hence comparing the resistance curves 
of these two cases essentially gives a method to extract 
the effect of the incompressible strips on the other layer. 

In order to investigate the relation between the in- 
compressible strips of the passive layer and the magneto- 
transport coefficients of the active layer qualitatively, we 



calculate the Hall resistances for a magnetic field interval, 
where QH plateaus are observed for both layers around 
filling factor two. In Fig.^Jwe show the Hall resistances 
(in units of the von Klitzing constant) vs. magnetic field 
for interacting (solid-lines) and non-interacting (dashed 
lines) systems. We start the discussion with the high- 
magnetic field regime (Q, c /Q,2 > 1.27, i.e right side of the 
vertical dotted line), which essentially corresponds to a 
density distribution similar to that shown in Fig. [SJu As 
there are no incompressible strips within both layers, no 
noticeable potential fluctuations are created due to local 
charge imbalance, therefore the Hall resistances of both 
layers are the same for interacting and non-interacting 
cases. If we examine the QH regime of the top layer (the 
regime between vertical dash-dotted line and dotted line 
in Figlllfl . it is seen that at least one incompressible strip 
is formed (cf. Fig. [SJd) , creating a potential fluctuation 
to the bottom layer. Meanwhile, the bottom layer is still 
compressible all over the sample, so the fluctuation can 
be screened nearly perfect, leading to a new distribution 
of the electrons and a small change in the Rjj curves, due 
rearrangement of the local conductivities. The situation 
is fairly different, in the interval 1.18 < Q c /fi2 < 1-25 
(the region between thin vertical solid-line and dashed- 
dotted line), where incompressible strips are formed at 
both of the layers (density distributions corresponding 
to Fig. IHt and Fig. EJi). In this regime, where both lay- 
ers produce fluctuations due to local charge imbalance 
and this fluctuation can not be screened perfectly ev- 
erywhere, which should result in a difference in the Hall 
resistance curves. This is observed for the bottom layer 
at the high-magnetic field edge, since the perturbation 
slightly shifts the maximum magnetic field value of the 
QH plateau. In this regime the quantized value of the 
Hall resistance of the top layer does not change, since 
it only depends on the presence of the incompressible 
strip. In the 1 < S1 C /S1 2 < 1.18 interval of FiglTTl there 
exists two stable incompressible strips in both layers sur- 
rounded by compressible regions, thus both layers are in 
the QH plateau. Here we should note that, although the 
potential fluctuations are created by the incompressible 
strips and are not screened perfectly, this perturbation 
changes only the positions or the widths of the incom- 
pressible strips, but not the value of the Rh within the 
plateau regime. Eventually it depends on the amplitude 
of the perturbation, i.e. if the amplitude is large enough 
to destroy the incompressible strips one does not observe 
the quantized value for Rh- The change in Rh can be 
observed for the top layer at somewhat smaller values of 
the magnetic field strength (0.8 < fi G /02 i$ 1-0) since 
the perturbation enlarges the plateau. In fact the fluc- 
tuation widens the incompressible strips, with respect to 
the non-interacting case and creates incompressible strips 
larger than the averaging length, \ av . Finally one ends 
with a wider plateau. The difference in the Hall resis- 
tance curves between the interacting and non-interacting 
systems, tends to disappear as the incompressible strips 
become narrower and move towards the edges by decreas- 
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ing the magnetic field. This is consistent with the pre- 
vious observation, that the amplitude of the fluctuation 
depends on the width of the incompressible region. 

In summary four magnetic field intervals are observed: 
(i) both layers are compressible and there exists no dif- 
ference for the Rh curves, calculated for the interact- 
ing and non-interacting case, (ii) top layer, at least, has 
an incompressible strip and creates a fluctuation, which 
shifts for the bottom layer the edge of the Hall plateau 
to higher magnetic field values, (iii) both layers show in- 
compressible regions, the perturbation generated by the 
bottom layer widens the incompressible strips of the top 
layer, leading to a wider plateau, (iv) the incompressible 
strips of both layers become narrower and move towards 
the edges and the fluctuation becomes inefficient, hence 
the difference is smeared out. It is useful to mention that 
the interacting system is an equilibrium solution, in the 
sense that the electrons of both layers rearrange their 
distribution until a full convergence is obtained (within a 
numerical accuracy), for a given external potential pro- 
file. 



V. COMPARISON TO THE EXPERIMENTS 

Here we report on the magneto-transport hysteresis 
observed in the bilayer systems measured at the Max- 
Planck Institut-Stuttgart by S. Kraus 2 -. In the experi- 
ments discussed here, the magneto-resistances are mea- 
sured as a function of the applied perpendicular magnetic 
field. The sweep direction dependence is investigated for 
matched and mismatch densities. Similar findings have 
already been presented in the literature 1718 , however, 
have been discussed in a different theoretical content. 

The samples are GaAs/AlGaAs double quantum well 
structures grown by molecular beam epitaxy (MBE). 
Separate Ohmic contacts to the two layers are realized 
by a selective depletion technique^. A technique devel- 
oped by Rubel et al.[35j was used to fabricate backgates. 
The metal gate on top of the sample acts as frontgate. 
The samples were processed into 80/im wide and 880^m 
long Hall bars. The high mobility samples are grown at 
the Walter-Schottky Institut. They have as grown den- 
sities in the range 1.5 — 2.5 x 10 11 cm -2 and mobility 
is 100 m 2 /Vs per layer. The barrier thickness is 12nm 
and the quantum wells are 15nm wide. The experiments 
were performed at low-temperatures (T 270mK) and the 
imposed current is always in the linear response regime 
(I~ 50nA). Other details of the experimental setup and 
samples could be found in an upcoming publication^. 

We show the measured Hall resistances of the top layer 
(solid lines for up and dashed lines for down sweep) and 
the bottom layer (dotted lines for up and dashed-dotted 
lines for down sweep) as a function of magnetic field 
strength and direction in Fig^J The data were taken 
at a sweep rate O.OlT/min and the base temperature is 
always kept at 270mK. Apparently it is seen that the re- 
sistances of the layers follow different traces depending on 



whether the sweep is in up or down direction, in certain 
B field intervals. We see that the sweep direction has no 
effect on the resistances at high magnetic fields (B > 8T) , 
where we would not expeclj^i to have any incompressible 
strips (see e.g. Fi^S^). The first exiting feature is ob- 
served when the bottom layer is in the plateau regime 
and the top layer is approaching to its plateau regime 
(6.5T< B < 8.5T). It is seen that the Hall resistance 
of the top layer follows two different traces depending 
on the sweeping direction, meanwhile the bottom layers 
resistance is insensitive and assumes the quantized value. 
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FIG. 12: The Hall resistances as a function of magnetic field 
for a GaAs electron bilayer system with a barrier width 12nm, 
considering two sweep directions. Data is taken at a magnetic 
sweep rate O.OlT/min. A strong hysteresis develops in Rh 
when the other layer is in a plateau. No significant effect is 
observed at low field regime. 

The plateau regime is extended up to B ~ 7.25T for 
down sweep and ends at B = 7T for up sweep, mean- 
while the bottom layer is also in the QH regime. It is 
easy to understand why the QH regime is wider for the 
down sweep considering the widths of the incompressible 
strips as follows: while coming form high magnetic field 
to low fields we first encounter with a wide incompress- 
ible region at the bottom layer (see for example Fig[5j3, of 
course here bottom layer has a higher density) that gener- 
ates a large potential fluctuation yielding a wide plateau. 
However, for the other sweep direction two narrow incom- 
pressible strips are observed first (Fig|8Ji) , which creates 
relatively weak fluctuations, that have a negligible effect 
on the plateau width. For smaller values of the magnetic 
field, we observe at the top layer that up and down sweep 
curves follow the same trace and have the quantized value 
at B = 6.8T. A similar hysteresis behavior is seen at the 
bottom layer resistance curves (5.7T< B < 6.5T), now 
the top layer is in the plateau regime for both sweeping 
directions. The resistances at low magnetic field inter- 
val (4T< B < 7T) do not show a prominent difference 
for two sweep directions. A less pronounced repetition of 
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the hysteresis is observed at the v = 2 plateau regimes, 
e.g. 3.4T< B < 3.6T and 2.9T< B < 3.2T. A common 
feature (also observed at different density mismatches) 
is that the hysteresis is seen at the active layer only if 
the passive layer is in the plateau regime. No hysteresis 
is observed if the QH regimes of the layers coincide (see 
discussion in Ref.|26|h We summarize the experimen- 
tal observations as follows: the hysteresis is observed in 
the active layer only if the passive layer is in the plateau 
regime. 

We should mentioned that the hysteresis effect is a non- 
equilibrium effect and " interacting" system is an equilib- 
rium solution, hence can not explain this effect directly. 
However, it is easy to attribute the widening of the top 
layers plateau (for down sweep) to the potential fluctua- 
tions created by the bottom layer, as already shown for 
the "interacting" case in Fig^] In accord with our nu- 
merical results (previous section) obtained for the "non- 
interacting" system, the measured quantities are inde- 
pendent of their history. We claim that the measured re- 
sistances are the analogues of the calculated resistances 
for the non-interacting system for the up sweep of the 
top layer, since the incompressible strips are narrow and 
the potential fluctuations are ineffective. 



ing potential. We compared the electron distributions for 
vanishing magnetic field and temperature with the finite 
ones, and observed that a local charge imbalance is cre- 
ated due to the formation of incompressible strips. These 
dipolar strips produce external potential fluctuations, as 
a function of applied magnetic field, to the other layer. 
We have investigated properties of this potential fluctu- 
ation by comparing the interacting and non-interacting 
systems for a few characteristic B values and obtained the 
Hall resistances by a local scheme proposed in our recent 
work«2i. We considered these fluctuations as a pertur- 
bation to the other layer and observed that they widen, 
stabilize and shift the plateaus as expected. In SecfVl 
we have report and attempted to obtain qualitative ar- 
guments for the hysteresis-like behavior by reconsidering 
the symmetry breaking of the sweep direction, based on 
the long relaxation times of the incompressible regions. 
This is done by considering the "interacting" and "non- 
interacting" schemes. Our results show that the Hall 
resistance curves follows different paths if both of the 
layers have incompressible strips within the sample. The 
amplitude of the deviation depends on the widths and 
the positions of the incompressible regions. 
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